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A Texture Atlas is an efficient color representation for 3D Paint Systems. 
The model to be textured is decomposed Into charts homeomorphic to 
discs, each chart is parameterized, and the unfolded charts are packed In 
texture space. Existing texture atlas methods for triangulated surfaces 
suffer from several limitations, requiring them to generate a large 
number of small charts with simple borders. The discontinuities between 
the charts cause artifacts, and make It difficult to paint large areas with 
regular patterns. In this paper, our main contribution is a new quasi- 
conformal parameterization method, based on a least-squares 
approximation of the Cauchy-Riemann equations. The so-defined 
objective function minimizes angle deformations, and we prove the 
following properties : the minimum is unique, independent of a similarity 
in texture space, independent of the resolution of the mesh and cannot 
generate triangle flips. The function is numerically well behaved and can 
therefore be very efficiently minimized. Our approach is robust, and can 
parameterize large charts with complex borders. We also introduce 
segmentation methods to decompose the model into charts with natural 
shapes, and a new packing algorithm to gather them in texture space. 
We demonstrate our approach applied to paint both scanned and modeled 
data sets. 
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Resources 



• An implementation is available in my OpenNL 
numerical library. 

• The 3D OpenSource modeler Blender now offers LSCM 
texture unwrapping. BrechtVan Lommel and Jens Ole 
Wund integrated OpenNL in Blender, more details 
given here and here . 

. Tutorials for using LSCM in blender: here , here and 
here (with video). 

• Microsoft's UVatlas (part of DirectX) uses our texture 
packing algorithm. 

• The 3D modeler Silo will use LSCM for texture atlas 
generation. 

• The 3D modeler Rheingold will use LSCM for texture 
atlas generation. 

Erratum and discussions 

In our paper about Least Squares Conformal Maps, 
presented at Siggraph in July 2002, we have erroneously 
claimed that the orientation of the triangles is preserved. 
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Further investigations have shown that the gradient of the 
LSCM function is equal to the gradient of the Dirichlet 
energy for inner vertices (i.e. Pinkall's cotangent formula). 
Therefore, like all other approaches based on this formula 
triangle flips may occur when the involved coefficients 
become negative. In practice, we have not observed this 
phenomenon in real data sets, and our segmentation 
method splits the corresponding neighborhoods. However, 
we have found the following 4-triangles surface, for which 
this problem occurs (Michael Floater has also a similar 
example on his web-page). 
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We are currently investigating ways of characterizing 
configurations where this problem occurs. Before we 
have a clean solution, this problem can be fixed using 
one of the following strategies: 

• As done by Desbrun and Alliez, use Rivara's method, 
enabling to remove obtuse angles by inserting new 
vertices In the triangulation. 

• As done by Sheffer, use constrained optimization to 
enforce the desired properties. In our case, 
constraining the sign of the Jacobian preserves the 
orientation of the triangles, 

. Make a chart boundary split the concerned 
neighborhoods. 

Note that this problem does not concern the other 
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properties (quadratic expression of the conformal 
energy, existence and uniqueness of the solution, 
independence to resolution, ability to extrapolate the 
border by constraining two vertices only, independence 
to a similarity applied to those two constrained 
vertices). 

Why did we make that mistake ? 

To study the properties of LSCM, we construct a mesh 
using two basic operations: 

• glue, i.e. add a new triangle to the border, by 

connecting three existing vertices; 
. join, i.e. add a new triangle and a new vertex to the 

border. 

Any triangulation homeomorphic to a disc can be 
constructed by these two operations. Then, we study how 
these operations modify the linear system solved by 
LSCM. For the existence and uniqueness (i.e. to prove 
that the matrix is of maximal rank), join is the easy case 
(adding a row to a matrix cannot decrease its rank) and 
glue requires more work (since it adds a triangle and a 
vertex, it adds a row and a column to the matrix). For the 
orientation of the triangles (i.e. the sign of the Jacobian), 
once we had solved what we had taken for the hard case 
(i.e. glue), we thought we were done. But for the 
orientation of the triangles, the easy and hard cases are 
reversed as compared to the rank problem !'. In fact, glue 
is the easy case, and join requires a more careful analysis. 
Knowing that, the proof in our paper shows in fact that a 
triangulation constructed by glue operations alone can be 
parameterized by LSCM without any triangle flip. 

DCP = LSCM 

Note: There is also a short note on Mathieu Desbrun's 
home page, indicating the equivalence between his 
Discrete Conformal Parameterization method 
(Eurographics, Sept 2002) and LSCM (Siggraph, July 
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2002). 

Mathieu's remark is based on the relation: Ec + A = Ed, 
where Ec corresponds to the conformal energy 
(minimized by LSCM), A to the area of the surface, and 
Ed to the Dirichlet energy. If A is fixed, then it is clear 
that minimizing Ec and Ed is equivalent. Since that 
{LSCM, DCP> minimize the {conformal, Dirichlet} 
energy of the function going from the 3D surface to the 
(u,v) parameter space, the area A under consideration 
is the area o f the (u.v) space_, therefore: 

. In the case of a fixed boundary, A is fixed, and 
DCP=LSCM (rem: this is no longer the case in 
presence of triangle flips) 

. In the case of a free boundary, A changes, and more 
precautions should be taken to determine whether the 
augmented DCP method (with Mathieu's additional 
terms to let the border free) is equivalent to LSCM. 
We have also shown the equivalence between both 
methods (the proof is too long to be included here). 
However, experimentally, we have obtained different 
results using the two methods, and did not manage to 
properly extrapolate the border of complex charts 
with DCP. We think this comes from a numerical 
conditioning problem: the LSCM formulation results in 
a symmetric matrix, which is not the case of DCP. In 
fact, LSCM and DCP are the same quadratic form, 
with the difference that LSCM is the Gauss-reducted 
expression of this quadratic form. 

Which method for which purpose ? 

• Fixed boundary on convex polygon: Floater's shape 
preserving parameterization and its derivatives, 
based on Tutte's barycentric mapping theorem 
(ensures the injectivity of the mapping, and is 
obtained as the resolution of a linear system). 

• Free boundary, small charts: Sheffer's Angle Based 
Flattening gives the best results, and preserves the 
orientations of the triangles. ABF also proposes 
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methods to solve global overlaps. The only limitation 
of ABF is speed (non-linear optimization problem) 

• Free boundary, large charts: I'd use LSCM ( guess 
why ;-) ), here are a few examples of large meshes 
parameterized using LSCM, The first image here is a 
single chart of 72,438 triangles !!. To my knowledge, 
LSCM is currently the only method capable of treating 
such a large chart with free boundaries. Note: we are 
currently developping ABF++ in cooperation with Alia 
Sheffer, an extension of ABF that can handle such 
large charts . 

• I like the results of the Stretch Optimization and the 
Signal Optimized Parameterization (Sander, Gortler & 
Hoppe), it is possible to use their optimization scheme 
as a post- processor for any parameterization method. 
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Abstract 

A Texture Atlas is an efficient color representation For 3D Paint Sys- 
tems. The model to be textured is decomposed into charts heme- 
omorphtc to discs, each chart is parameterized, and the unfolded 
charts are packed in texture space. Existing texture atlas methods 
for triangulated surfaces suffer from several limitations, requiring 
them to generate a large number of small charts with Simple bor- 
ders. The discontinuities between the charts cause artifacts, and 
make it difficult to paint large areas with regular patterns. 

In this paper, our main contribution is a new quasi-con formal pa- 
rameterization method, based on a least-squares approximation of 
the Cauchy-Riemann equations. The so-defined objective function 
minimizes angle deformations, and we prove the following proper- 
ties: the minimum is unique, independent of a similarity in texture 
space, independent of the resolution of the mesh and cannot gener- 
ate triangle flips. The function is numerically well behaved and can 
therefore be very efficiently minimized. Our approach is robust, 
and can parameterize large charts with complex borders. 

We also introduce segmentation methods to decompose the 
model into charts with natural shapes, and a new packing algorithm 
to gather them in texture space. We demonstrate our approach ap- 
plied to paint both scanned and modeled data seta. 
CR Categories: 1.3.3 [Computer Graphics) Picture/Tmagc Gen- 
eration; L3.5 [Computer Graphics]: Three-Dimensional Graphics 
and Realism—Color, shading, shadowing and texture; T.4.3 pmoge 
processing]: Enhancement — Geometric Correction, Texture 
Keyword A; Texture Mapping. Paint Systems, Polygonal Modeling 
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1 INTRODUCTION 

A 3D paint system makes it possible to enhance the visual nppoa *- 
ance of a 3D model by interactively adding details to it (color i, 
bump maps . . .). If the discretization of the surface is fine enough, 
it is possible to directly paint its vertices [1]. However, in mo it 
cases, the desired precision for the colors is finer than the geome:- 
ric details of die model. Assuming that the surface to be painted 9 
provided with a parameterization, it is possible to use texture ma] >- 
ping to store colors in parameter space {9). Parametric eurfocis 
(such as NURBS) have a natural parameterization. For other jCpp> 
se motions, such as polygonal surfaces, rinding a parameterizatk n 
is non-trivial. To decorate polygonal models with regular patterns, 
the lapped textures approach [26] can be applied: local overlappirg 
parameterizations arc used to repeatedly map a small rex rune swatt h 
onto a model, 

A texture atlas is a more general representation (sec, e.( 
[13, 20, 23]). The model to be textured is partitioned into a s?t 
of pares ho meom orphic to discs, referred to as charts* and ca< h 
of them is provided with a parameterization. A texture atlas ci n 
be easily represented by standard file formats and displayed usii g 
standard texture mapping hardware. When used in a 3D paw* sys- 
tem, a texture atlas should meet the following requirements: 
• the chart boundaries should be chosen to minimize texture ani- 

fects, 

« the sampling of texture space should be as uniform as possible, 
■ the atlas should make an optimal use of texture space. 

The generation of a texrure atlas can be decomposed into the 
following steps: 

1. Segmentation: The model is partitioned into a set of charts. 

2. Parameterization: Each chart is 'unfolded', i.e. put in con> 
spondence with a subset of 3t 2 . 

3. Packing; The charts arc gathered in texture space. 

The remainder of this section presonts the existing methods lar 
these three steps, and their limitations with respect to the requue- 
mcnts mentioned above. We then introduce a new texture atlas gen- 
eration method, meeting these requirements by creating charts wtfi 
natural shapes, thus reducing texture artifacts. 
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1.1 Previous Work 

Segmentation into charts. In [14] and [23], the model is interac- 
tively partitioned by the user. Tb perform automatic segmentation, 
Maillot et al. [201 group the facets by their normals. Several multi- 
resolution methods [7, 16] decompose the model into charts corre- 
sponding to the simplicca of the base complex. In [27], Sander et 
al "use a region-growing approach to segmentation, merging charts 
according to both planarity and compeotness criteria. All these ap- 
proaches are designed to produce charts that can be treated by ex- 
isting parameterization methods, which arc limited to charts; with 
convex borders. For mis reason, a large number of charts is gen- 
erated, which introduces marry discontinuities when constructing a 
rexturc alias. 

Chart parameterization. Discrete Harmonic Map, described by 
Eck et al [4], are the most widely used. They arc approximations 
of Continuous Harmonic Maps [5], minimizing a metric dispersion 
criterion. Pinkal and Polthier [24] have shown the link between this 
criterion and another one named conformality, and have expressed 
both in terms ofOirichht energy. Haker */ al [8] describe a similar 
method in the specific case of a surface triongulation homeomor- 
phic to a sphere. 

The theory on graph embedding has been studied by Tuttc [30], 
where Barycentric Maps are introduced- The bijectivity of the SO- 
defined parameterization is mathematically guaranteed. Floater 16] 
proposes specific weights improving the quality of the mapping, in 
terms of area defbrmations and conformulity. In [18], a method is 
proposed to take additional constraints into account. 

In all the methods mentioned above, since conformality is ex- 
pressed as an indirect coupling between the parameters, boundary 
conditions arc required, i.e. boundary nodes need to be fixed on 
a convex border in parameter space. Other expressions of confor- 
mality. such as the non-linear MIPS method [10], make it possi- 
ble to overcome this problem, and let the boundary nodes be free 
to move, However, this latter method requires a time-consuming 
non-linear optimization, and may get Stuck in a local minima of 
the non-linear function. In [12], Hurdal et al propose a method 
based on circle packings, which are cortain configurations of cir- 
cles with specified pattern of tangencies known to provide a way 
to approximate a conformal mapping. Building circle packings is 
however quite expensive. The approach proposed in [28] consists 
in solving for the angles in parameter space. Tt result* in a highly 
constrained optimization problem. Other methods [17, 25, 31] can 
also extrapolate the border, but do not guarantee the absence of tri- 
angle flips and require interaction with the user. We introduce here 
a conformal mapping method, offering more guarantees, efficiency 
and robustness than those approaches. 

In the case of texture mapping, not only the bijectivity of the 
parameterization should be ensured, but also its ability to make an 
optimum use of texture memory, and to accurately represent a sig- 
nal stored in texture space. Sander et. al. [27] describe an approach 
to minimize both a texture stretch criterion, and texture deviation 
between level of details. Since their approach is independent from 
the initial parameterization method, it can be applied to optimize 
the sampling of the parametcrizations constructed by our method. 

Charts packing Id texture space. Finding the optimal packing 
of the charts in texture space is known ae the bin packing problem. 
It has been studied by several authors, such as Milenkovic (see, 
e.g., [21 ]), but the resulting algorithms take a huge amount of time 
since the problem is NP-complcte. To speed up these computa- 
tions, several heuristics have been proposed in the computer graph- 
ics community. In the case of individual triangles, eueh a method in 
described by several authors (see, e.g., [3]). In the general case of 
charts, Sander et al [27] propose an approach to pack the minimal 



area bounding rectangles of die charts. In our case, since the chart 
can have arbitrarily shaped borders, the bounding rectangle can b< 
tar away from the boundary of the charts. Therefore, a lot of texrun 
space can be wasted. For this reason, we propose a more accuratx 
packing algorithm that can handle Che complex charts created b} 
our segmentation and parameterization methods. 

1.2 Overview 

The paper is organized OS follows. Since it is our mam contribu- 
tion, we will start by introducing Least Squares Conformal Map i 
(LSCMs), a new optimization -based parameterization method with 
the following properties (see Section 2 and Figure I): 
« Our criterion minimizes angle deformations and non-uniform 
scaling*. It can be efficiently minimized by classical MA meth- 
ods, and does not require a complex algorithm such as the one i 
used in [12] and in [28]. 

• We prove the existence and uniqueness of the minimum of 
this criterion. Therefore, the solver cannot get stuck in a locr 1 
minimum, in contrast with non -linear methods [10, 25, 27, 31] 
where this property is not guaranteed. 

m The borders of the charts do not need to be fixed, as with mojt 
of the existing method* [4. 6, 18]. Therefore, large charts witi 
arbitrarily Shaped borders can be parameterized- 

• Wc prove that the orientation of the triangles is preserved, whic l 
means that no triangle flip can occur. However, as in [28], over- 
laps may appear, when the boundary of the surface sclf-intersecis 
in lexture space. Such configurations are automatically detected 
and the concerned charts are subdivided. This problem was se - 
dom encountered in our experiments (note that as with classic:. I 
methods [4, 6], if all the border nodes are fixed on a convex pol) - 
gon. no overlap can occur). 

m We prove that the result is independent of the resolution of tbe 
mesh. This type of property may be usefull to reduce tcXTUie 
deviation when parameterizing different level of details of tie 
same object. 

In Section 3, we present a new segmentation method to decom- 
pose the model into charts. Thanks to the additional flexibility of- 
fered by LSCMs, it is possible to create large charts correspond- 
ing to meaningful geometric entities, such as biological features nf 
characters andanimals. The required numbor of charts is dramaii- 
cally reduced, together with the artifacts caused by the discontinu- 
ities between the charisi. Moreover, thcue large charts facilitate the 
use of regular patterns in a 3D paint system. 

Section 4 presents our method to pack the charts in texture spac i. 
Since our segmentation method can create charts with complex her- 
ders, we pack the charts more accurately than with bounding ret t- 
anglcs, as in previous approaches. Our method is inspired by the 
strategy used by a 'Tctris* player. The paper concludes with some 
results, on both scanned and modeled meshes. 

2 LEAST SQUARES CONFORMAL MAPS 

In this section, wo focus on The problem of parameterizing a chi rt 
homeomorpbic to a disc. It will then be shown bow to decompose 
the model into a set of charts, and how to pack these charts in texture 
space. 

2.1 Notations 

• scalar* arc denoted by normal characters x, y, w, v\ 

• vectors are denoted by bold characters x = (x> y). 

• complex numbers arc denoted by capitals U = (u + iv). 

• vectors of complex numbers arc denoted by bold capitals U, 

• maps and matrices are denoted by eursivo fonts W, X. 
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Figure I : The body of the scanned horse is a text case fir the robustness of fa method; it to a tingle very largo chart <?T a confite bor*r. ******* 

iso~par™*ter curves; B: The corresponding unjoin surface, where the bonier has ben *ut*»otlcaify extrapolated: C: Tftw- am ^ lAr xv^« e^lwfai* to a db*; A. TV. r 
parumderttation is rvbvtt. and not qffecrcd by the targe, triantfes in the circled area (caused by shadow zones appearing during the scanning process). 



2J2 Conformal Maps 

Tn this section, we quickly introduce the notion of conformal map. 
We will present further a new way to approximate the conformality 
criterion and the mathematical properties of thia approximation. 




Figure 2: tn a conformal map. the torment wtorx to the hn-u and w the L\a*v curves 
are orthogonal and have the same Ungth. 

As shown in Figure 2, an application X mapping a (at, v) domain 
to a surface is said to be conrormal if for each (u,v), the tangent 
vectors to the iso-n and iso-v curves passing through ^(u, v) are 
orthogonal and have the same norm, which can be written as; 



dx f . 



*7 



0) 



where N(u, v) denotes the unit normal to the surface. In Other 
words, a conformal map is locally isptrvpic, i.e. maps an elemen- 
tary circle of the (u, v) domain to an elementary circle of the sur- 
nice. 

It is possible to rewrite Equation 1 using differential operators, 
such as Laplace-Beltrami.asdonein [24] and in [8], which results 
in the well known cotangent weighting coefficients (sec e.g. [4]). 
The (u, v) parameters are then found to be the solution of two sep- 
arate linoar systems, one for i* and ono for v. The relation between 
u and v is indirectly taken into account by the right hand sides of 
the two systems. For this reason, this type of method requires the 
border to be fixed on a convex polygon. The MIPS method [10] 
does not have this restriction, and expresses conformality as a re- 
lation linking the coefficients or the metric tensor. However, the 
resulting equations are non- linear. Another approach has been de- 
scribed In [28], based on Ihc remark thai the criterion defining a 
conformal mapping should be independent of a translation, rotation 
and scaling in parameter space (i.e. a similarity). The unknowns 
are the angles at the corners of the triangles. This requires a time- 
consuming constrained optimization method. 

Rather than discretizing the Laplace operator at the vertices of 
Che tri angulation, we instead take the dual path of considering the 
conformaJity condition on the triangles of the surface. Using the 



fact that a similarity can be represented by the product of compter 
numbers, we show how to turn the conformality problem into a * 
unconstrained quadratic minimization problem. The u and v pi - 
ramctcrs are linked by a single global equation. This direct coi - 
pling of the u and v parameters makes it possible to efficiently pi - 
rameierize largo charts with complex borders, as ahown in Figuie 
1 . Tn this example, the cuts have been done manually (Figure 1-C \, 
to create a large tctit cose for the robustness of the method, (It will 
be shown in Section 3 how to automatically cut a model into char s 
homeomorphic to discs.) 

RiemaiuVs theorem states that for any surface S homeomorpn c 
to a disc it is possible to find a parameterization of the surface sa 
isfying Equation 1. However, since we want to use the resulting 
parameterization for texture mapping, we add the constraint th it 
the edges of the triangulation should be mapped to straight line 3, 
and the mapping should vary linearly in each triangle. With th s 
additional constraint, it is not aJways possible to satisfy the cortfo*- 
mality condition. For this reason, we will minimize the violation uf 
Ricmann's condition in the least squares sense, 

2.3 Conformality in a Triangulation 

Consider now a triangulation Q = {[1 . . ■ n], T> (Pj)i*>«n K 
where [l . . .n],7» > 3, corresponds to the vertices, where T is 
a set of ri triangles represented by triples of vertices, and whe e 
Pj e denotes the geometric location at the vertex j. We su i- 
pose that each triangle is provided with a local Otthonormal basis, 
where (a-i.s/i), ya) arc the coordinates of its vertic *s 

in this basis (i.e., the normal is along the z-axis). The local bas ss 
of two triangles sharing an edge are consistently oriented. 

We now consider the restriction of X to a triangle T and apply 
the conformality criterion to the inverse map U : (», y) *-* (v, 0 
(i.e. the coordinates of the points are given and we want their pa- 
rameterization). In the local frame of the triangle, Equation I be- 
comes 

dX OX ft 

where X has been written using complex numbers, i.e. X = x + t y. 
By the theorem on the derivatives Of inverse functions, this impli ts 
that 

<*) 



where U = u + iv. (This is a concise formulation of the Cauchy- 
Riemann equations.) 
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Since this equation cannot in general be strictly enforced, we 
minimize the violation of the confonnalitv condition in the least 
squares sense, which defines the criterion C: 



&i\ 2 , A law awl 3 . 



where At is the area of the triangle and the notation |z| stands for 
the modulus of the complex number z. 

Summing over the whole tri angulation, the criterion to minimize 
is then 

rcr 

2.4 Gradient in a Triangle 

Our goal is now to associate with each vertex j a complex num- 
ber Uj such thai the Cauchy-Rtemann equation is satisfied (in the 
least squares sense) in each triangle. To this aim, let us rewrite the 
criterion C(T), assuming the mapping U varies linearly in T. 

We considor a triangle {(x i ,y\) r («a , ya), (*3, V*)} of R 2 . with 
scalars U\ , , associated with its vertices. We have; 

where rfx - O^I/a - $1*2) + (*22/s - 1/2*3) + (rcayi - itjii) is 
twice the area of the triangle. 

The two components of the gradient can be gathered in a com- 
plex number: 



|a + = ± rift w* iv 3 )(ui U2 m) r , 

ox oy flr 



where 



The Cauchy-Riemann equation (Equation 2) can be rewritten as 
follows: 

t/y OtT* 

where ^ - iij + *Vj. 

The objective function thus reduces to 

C{V = (Ui U») T )= E C{T) % with 

where triongle T has vertices indexed by j\, ja, jV (We have mul- 
tiplied C(T) by a factor of 2 to simplify the expression.) 

2.5 Least Squares Conformal Maps 

C(U) is quadratic in the complex numbers U\ C/ n . so can be 

written dawn as C(U) = U'CU, 0> 

where C is a Hermitian symmetric ra x n matrix and the notation 
IT stands for the Horminan (complex) conjugate of TJ. C is an 
instance of a Hermitian Oram matrix, i.e. it can be written as 

where M « (my) is me sparse n' x n matrix (rows arc indexed 
by triangles, columns arc indexed by vertices) whose coefficient is 



if vertex j belongs to triangle Ti, 
otherwise. 




FigUTC 3: Our LSCM parameterization is insensitive resolution of the men) . 
lit tso-fiarojuotur tvrtttx obtained on «i waniv modi (W\ptn A) and tin a fine tmt 
(Figure B) an Identical, and remain stable when the resolution varies within a mes 1 
(circled zone In fttiures C and O). 

For the optimization problem to have a non-trivial solution, som 5 
of the UiS must be set to a priori values. Let us decompose the vec - 
tortJas (U/,Uj) T , where U/ is the vector of/nec coordinates <f 
TJ (the variables of the optimization problem) and V F is the vectcr 
of pinned coordinates of U, of length p (p ^ n). Along the sams 
lines, M can be decomposed in block matrices as 

M - (M f M p ) , 

where Mf is n ri x (n - p) matrix and M P is a ri x p matrh.. 
Now, Equation 3 can be rewritten as 

<7(U) = V m M'MTJ = WMUf = \\M/Vf + M P V P \\ 2 , 

where the notation ||v|| 2 stands for the inner product < v, v > (7 
stands for the conjugate of v). 

Rewriting the objective function with only real matrices and ve» - 
tors yields 



C(x) = ||-4x-b|| 3 ( 



(<•) 



with 



*-W (mi -%)(%)■ 

where the superscripts 1 and 2 stand respectively for the real ai d 
imaginary part, ||v|| stands this time for lie traditional Z^-norm »>f 
a vector with real coordinates and x = (U} T ,U} )~ is the vect >r 
of unknowns. a / 

Note that A is a 2rt' x 2(n - p) matrix, b is a vector of R n ai d 
x is a vector of H 2(r1 " p) (the im and Vi coordinates of the vcrtic js 
in parameter space that are allowed co move freely). 

2-6 Properties 

Tho above minimization problem has several fundamental proper- 
ties which are proved in the appendix: 

• The matrix A has full rank when the number of pinned vertici s, 
i.e. p, is larger than or equal to 2. 

• As a consequence, the minimization problem has a unique sol i- 
tion when p > % given by x = (A T A)~ l A T h, The best vol ic 
for p is 2, since in this case the mapping hi can be fully confonral 
if the surface is developable (i.e. the minimum of the objecti fc 
function is zero). In our experiments, we have pinned the t\ k> 
vertices maximizing the length of the shorted path "between them 
(i.e. the graph diameter). 

• The solution to *hc minimization problem is invariant by a simi- 
larity in texture space. 

• Tho solution to the minimization problem ia independent of I \c 
resolution of the mesh. This property is illustrated in Figure 3 

• In texture space, all the triangles are consistently oriented if uc 
pinned vertices are chosen on the boundary of T. In other won Is, 
triangle flips cannot occur. 
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Figure 4: A. Raxult of thr footurr* detection algorithm; B: The di*t aT u*_to_sced Junction it not an optimal choice for driving our chart growing prvee**; C: Thi 
distort jojeatur** function shows ho-tontounwlth mo* natural thapes: to it**»tt of aur segmentation algorithm, *> the dUtawcJoJtztturctpnvdpn, 



3 SEGMENTATION 

The segmentation algorithm decomposes the model into a set of 
charts. The design of the algorithm aims £it meeting the following 
two requirements: 

1. charts boundaries should be positioned in such a way that most 
of me discontinuities between the charts will be located in zones 
where they will not cause texture artifacts, 

2. charts must be homeomorphic to discs, and it must be possible to 
parameterize them without introducing too much deformation. 

For the first point, since the shading models depend on the nor- 
mal, zones of high curvatures cause sharp variations of lighting. In 
these zones. a texture artifact will not be noticeable, since it will 
be negligible compared ro the shading variation. Thus, to minimize 
artifacts, wo will design the segmentation algorithm in such a way 
as to avoid chart boundaries in flal zones. In other words, it is suit- 
able to generate large charts with most of their boundaries in high 
curvature zones. 

For the second point, we will present an automatic approach, 
mimicking the way a user manually sogments a model. Basically, 
the user attempts to decompose the model into parts resembling 
cylinders. Tb detect such cylinders, we will use an approach in- 
spired by Morse theory, characterizing a function defined over the 
surface (see, e.g , [29]). 

The next two sections present a feature detection algorithm, 
which finds curves corresponding to high eurvaiure zones of Ihc 
model, and a chart growing algorithm making them meet at these 
feature curves. Then, it will be shown how to validate the result, 
and subdivide the charts if needed. In the remainder of this section, 
we suppose that the surface is represented by a halfedgc based data 
structure (see, e.g., [19]). 

3.1 Detect Features 

The features detection phase can be outlined as follows: 
I Compute a sharpness criterion on the edges. We use here the 
second order differences (SOD), i.e. tho angle between the nor- 
mals, as in [1 1 ]- It is also possible to use more elaborate criteria. 
2. Choose a threshold r SO that a certain proportion of the edges is 
filtered out In our examples, we kept 5 percent of the detected 



3. For each of me remaining edges, grow a feature curve by apply- 
ing Algorithm 1 . 

Algorithm 1 attempts to anticipate the best paths, and filters out 
ihc small features caused by noise. Tagging the neighborhooda of 
the detected features avoids generating a large number of features 
in zones of high curvature. In our examples, the parameters arc set 



expand_reat«rc_curve<halfedgc atari) 

for hnlfedge h 6 { a tart, opposite* **<nri> } 
halfedgc fi' — h 
dn 

use dcpMlKi search to find the siring S of hattedges 
Hkirtms Wiih ttrtd * uch *ant: 

• two consecutive ho Hedges of 5 share a vertex 

« ihc length of S < dian max_atri7>aj«nath 

• 6harpnex»(S) — r± e s aharpnttm^) Is maximum 

• no halfedgc of S goci backward (relative to V) 

• no haU'edge of S Is tagged v5 u fbniuro neighbor 
h' — second hem of 5 

append />' to dr.tr.r.Ltid^anturG 
whil*(tiha.rpnBBB(S) > m.*x_atringjcngtk X t) 
end //for 

lt(langth{dctccted_Jtin.*u.r&)> min, _J catvrr.jMgth) rticn 
tag ihe elemenui of detcct*d_f mature as feoturw 
tag ihc haJredgeM in ihc neighborhood of dc.tcc* nd^f nature 
is tenure neighbors 

tml/S expand Jeatimjcvn* 

Algorithm 1: Features growing 



as follows: max string length = 5, which controls the size »»f 
the discontinuities to be fiTlcd, and rnin _JeatureJieJigth ~ 15. 

3.2 Expand Charts 

Once the sharp features have been detected, the charts can be cr > 
ated. Our method is a greedy algorithm, expanding all the charts 
simultaneously from a set of seeds. It is similar to the «^soun« 
Dijksrra algorithm used in [4] and to the region-growing paradig n 
used in computer vision. Since we want chart boundaries* to meet it 
the level of features, the a-source algorithm is modified as follows: 

• To select the set of seeds, the intuitive idea it* to 'reverse engi- 
neer* the expected result. More precisely, we use the follow- 
ing method: a front i# propagated from the borders and the fea- 
ture curves detected by the previous algorithm, to compute a 
distance_to_jcatures function at each facet Then, the seels 
are found to be the local maxima of this <&Btance_to^f eatun :b 
function, 

• For closed surfaces without any detected feature, propagation is 
initialized from the two extremities of a diameter of the facets 
graph, as done in (15]. 

• Our s-sourcc propagation uses — distance JLo^fev-twe* as t'ic 
priority function, rather than di$tanceJ,o_aeeds. The adva i- 
tagc of this approach is shown in Figure 4. 
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Figure 5: A; Out segmentation algorhhm defect* cylindrical vhap**: B< An nddh 
lioruil cut is added to 'sock-shaped' extremal cytinders. 



cxpH ri <l_ch arts 

P Hariiy_qucuo<Juilfodgb> Nr.np jwrtdd by rlint(ftkpnr.(hal/r.ttnc)) 
on<cd&> cKort_bov.nda.rieB initialized with nil ih* «dgea of ihc surface 
// Initialize Heap 

fa reach fee*i F where d 1 51(F) is a locul maximum 

create a new chan wixh seed F 

add this hnlfedfi« of F k> H top 
Bad //forrach 

// Charfj-srv^ins phase 
whMc{//Kop is not cmpiy) 

hulfedso h — c £ Heap such lh«l dlsi(e) is nuinlmum 

remove h from J/eop 

facer F — ftccl(/0 

fcac F„ pp «- ihc opposite facet of F relative to h 
|f( ehnr^F^pp) is undefined ) then 

add F OJip to ch*n(F) 

remove F flnom chart Jmundarixis 

renMve non-cxircmnl odgoi from ci\artji6-ndarii:* t 

//(I.e. e«fcw thai do not link t»* ttther chart boundary edges) 

odd the hairodgca of F„ f , T belonging to 
chart JboundariAa to Heap 
Clioir(chan(F np "|,)^chafi(F) nad 

mait_d!ai(chttrt(F^) - dlst(F) < c and 

max"dist(ehart(F„ r p» - disuTF) < e ) then 

merge chwtfF) and eiutrtfFop,,) 
enA//lf 
na//whlh 
end /> apandjpharts 

Algorithm 2: Charts growing. 



• Charts are merged if they meet at a small distance 
d < e from their Good. Tn our ex peri merits, e = 
riXO-xdist/4*, Where maxdiat denotes the global maximum of 
distajicejto _f matures. 

Our charts growing algorithm (Algorithm 2) uses the following 
data: 

• distance Jo_featwre Ls stored in each facet F. and denoted 
dist{F)\~ 

• for each chart C, the scalar maxjiiat{C) denotes the maximum 
distance to features for all the facets of C\ 

• The set of edges chart J><nrndarie9 represents the borders of 
all charts. It makes it possible for a patch to be its own neighbor 
while remaining a topological disc. 

As shown in Figure 5, our algorithm can detect cylindrical 
shapes, as the approach proposed in [15]. In our case, two con- 
figurations can be distinguished: 

• The cylinder corresponds to on extremity, such as the 'ringers* 
of the dinosaur's wings (Figure 5-B). The corresponding charts 
have the shape of a sock. This configuration is detected by corn- 




Figure 6: A: The dinosaur's head made a/a single chart; «V Despite dw absence < f 
Iripntfe flips, overlaps may occur, by self-intcmxtbms of the border, O The t 

con be ivmoved by subdividing die chart. 



miting the area/perimeter ratio. In this case, to facilitate the pa- 
rameterization phase, a cut is added by starting from the seed an i 
cutting the edges along the steepest-descent path. 
• The cylinder is non-oxtremal, and therefore non-capped. Algc - 
rithm 2 generates suitable boundaries, without requiring any spe- 
cial treatment (the resulting chart is 'rolled' around the cylinder ■. 



3.3 Validate Charts 

The so-constructed charts art then parameterized using the method 
presented in Section 2. After that, the following two criteria aie 
tested: 

• As mentioned in Section 2. no triangle flip can occur, and Ire 
border can be extrapolated. However, since the border can re 
non-convex, a new class of overlaps can be encountered. The y 
are caused by a self-intersection of the border* as In [28]. Sue h 
configurations can be efficiently detected by the hardware, t y 
drawing the parameter space in stencil mode. The stencil pixe-s 
drawn more chan once correspond to overlaps. If such ovcrlai is 
are detected, the corresponding chart is subdivided, by cutting it 
along the edges on the border of the overlapped zone, as shov. n 
in Figure 6. Note that our segmentation algorithm would n it 
generate a single chart for the dinosaur's head (see Figure 5). n 
practice, the overlap problem has seldom appeared in our expe r- 
imems. and was caused by tiny loops Formed by the border. 

• Our criterion respects angles very well, as shown in the results 
section. As far as areas are concerned, large charts with zones • >f 
high curvature may result in large area variations, lb detect the le 
problems, the minimum and maximum model area/texture ar» 
ratio is measured over the facets, If the max/min ratio is great sr 
than a certain threshold, the concerned charr is split, by growii ig 
two charts from the facets corresponding to the minimum and the 
maximum. In the examples shown here, the threshold has bom 
set to 2. 

4 PACKING 

Once the model is decomposed into a set of parameterized charts, it 
is possible to create a texture atlas by merging all the (u, v) domai is 
of the charts. Usually, only a limited amount of texture memory is 
available. It is then suitable to minimize the unused space. In other 
words, given a set of possibly non-convex polygons, we want to 
find a non-overlapping placement of the polygons in such a ww 
that (he enclosing rectangle is of minimum area. The so-obtain ;d 
texture coordinates are then rc-scalod to fit the size of the texture 

The packing problem is known to be NP-complctc (see, e.g„ [2 1 ] 
and [22]). Approaches based on computational geometry ahiw 
good performances in terms of minimization of lost area, but trc 
not efficient enough for large and complex data sets. For this rea- 
son, several heuristics have been proposed in computer graphbis. 
For instance, in the method proposed by Sander et at. [27], tae 
bounding rectangles of the charts are packed. In our case, since 
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Figure 7: A : Ow packing algorithm Insert the charts one by one, and maintains the 
'horizon ' (In blue) during the prmv- Each chart (In preen) Is inserted at the position 
mtnimtslng the 'wasted space' (In Mack} between hs bottom horizon (In pink) and tiie 
current horizon. The top horizon (in red) of the current chart b ihtrt used ro vpdm the 
harfian. 8: Result on the dtaoxavr data sot. 

the border may have an arbitrary shape, the bounding rectangle is 
not an accurate approximation. 

For this reason, we propose a different algorithm, that packs the 
charts directly rather than their bounding rectangles. As in [2], our 
algorithm is inspired by how a "Tetris' player would operate, but 
without approximating the charts by their bounding boxes: 

1 . Each chart is repealed to make its area in (u, v) space equal to 
its area in (x, y, z) space. 

2. The maximum diameter of the charts are oriented vertically and 
sorted in decreasing order. 

3. As shown in Figure 7. for each chart C» the top horizon hj. (in 
red) and bottom horizon h£ (in pink) is computed. Rather than 
being represented by their bounding rectangles, the charts are 
approximated by the area between the two corves tic and h&. 
As in classical approaches, in order to avoid unwanted blends 
caused by mip-mapping, an additional margin is added to the 
horizons. 

4. The charts ore inserted One by ©nc> using the method described 
below. 

As shown in Figure 7-A, the piecewise linear funciion h(u) rep- 
resenting the 'horizon 1 is maintained by the algorithm. For each 
chart C, the uc coordinate at the lower left corner of C is cho- 
sen in such a way that the lost space (in black) between the bottom 
horizon h& of C (in pink) and the current horizon h (in blue) is 
minimized (see Figure 7-A). Then, *e horizon h is updated using 
the top horizon he of the current chart (in red). Since the parame- 
ter space is discretized into texels, it seems natural to represent the 
horizons by arrays of discrctized values, with texture resolution. 
This makes the algorithm much simpler than using piecewise linear 
functions. For a chart C, all discrete values for uc are tested. The 
algorithm performs well, and takes less than one second to process 
all the data sets we have tested. 



5 RESULTS 

Wc have applied our method to different data sets, comprising 
meshes created with o 3D modeler (using subdivision surfaces) and 
scanned meshes. The results arc equivalent to those obtained with 
MIPS [10], but with mathematical guarantees, and can be more ef- 
ficiently computed. As shown in 'Table I, the stretch (sec [27]) 
measured on tbe result is of the same order as when using stan- 
dard methods <e.g. [4, 6]). To optimize the mapping, it is easy to 
post-process the result of our method by the algorithm proposed in 
p7]. Since the border nodes are natural ry positformed (rather than 
arbitrarily fixed on a convex polygon), the result can be better (see 

^Figure 8 shows some texture atlases. Note theprcsence of small 
charts, most of them corresponding to geometric details of the mod- 
els (teeth, hoofs . . .). This is not a problem for most paint systems, 



' IP LAW ■> PTO 121016/019 





Harmonic Maps 




stretch (before optim.) 


3.2 


3.5 


strctcn (arter optim.) 


1.65 


L52" 



Table 1 1 Stretch optimi=a(fart afttic 'cow head ' data net 




Figure 8 : Aotu /teii and associated tabu* space constructed by out method. 





dinosaur 


"skull 


Bunny 


none 


tt vertices 




~To7949" 


34334" 


4feVS5- 


d facets 


14JW4 


15,124 


&W51 


96/. 68 


1J charts 


43 


40' 


23 


44 


segmentation time (s; 


H 


17 


30 


43 


parameterization time (sj 


11) 


"23 




190 


packing rano (rectangiesj 


0.4K 


'0.51 


0.43 


■0.3' 


packing ratio (our algo.) 


"0.55 


"035 


0.6 


0.51 


stretcn (oerore optim.) 


-2.9 


2.5 


Kid 


U'. 


sire ten (after optim.) 


1/2* 


1.55 " 


-1.14 





Table 2 : Statistics and timings. 
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Figure 1 0: Hand-painted 3D models. Our !,$CM method fi>cilH°m the use of procedural textures and compt** patterns. 




Figure 9: Angle and ana deformmtom histograms ('Horse ' data setf. 



that can treat them property. Some examples of textured models 
are shown in Figure 10. Table 2 shows the sizes of the dam sets, 
the number of created charts, and the following statistics, obtained 
on a 1.3 GHz Pentium III (note thai the timings for the packing 
algorithm are not included, since they are negligible): 

• time to segment the model into charts', 

• time to parameterize the charts. Our LSCM criterion (Equation 
4) is minimized using me CG (Conjugate Gradient) algorithm. 
The independence to resolution suggests that a multi-grid ap- 
proach would be even more efficient; 

• packing ratio obtained using an enclosing rectangle packing ap- 
proach [27] and our algorithm. 

- stretch measured before and after applying Sander el. at.'s opti- 
mization method as a post-processing (see (271). 

The left histogram in Figure 9 shows the distribution of the an- 
gles in degrees between u and v gradient vectors. The mapping is 
nearly conformal in each triangle (the differences of lengthes be- 
tween the ti and v gradients we have measured arc very near to 
zero). The right histogram shows the area deformations obtained 
with the 'Horse' data set, before stretch optimization. This his- 
togram, showing texture area/model area ratios has been normal- 
ized, i.e. scaled in such a way that the mean value is mapped to 
I . Note that since the mapping is nearly isotropic in each triangle, 
the I? and L°° stretch histogram (not shown here) have exactly 
the same appearance as the area histogram. As can be seen, even 
though our LSCM criterion is not designed to punish area defor- 
mations, few facets are distorted, and can easily be fixed by post- 
processing using Sander et. af.'s method The resulting texture at- 
lases combine the advantages of LSCM (few chart discontinuities) 
and stretch-optimized parameterization (uniform sampling). 



CONCLUSION 

In this paper, we have presented a new automatic texture atlas gei - 
eration method for polygonal models. Overall, we have proposed a 
complete and mathematically valid solution to the parameterization 
of complex models which proved to he more efficient and robu-il 
than existing methods and available tools in real production env - 
ronmcnts. Our segmentation algorithm, driven by detected featurt s 
and inspired by Morse theory, decomposes the model into chars 
vrith natural shapes, corresponding to meaningful geometric ent - 
tics. These two algorithms may have applications in otheT domain 
such as re- meshing and data compression. We have successfully 
applied our technique to both scanned and synthetic data sees, mail- 
ing it possible to use existing 3D paint systems with them ()0ee] * 
Pmnt3D, Painter). In future works, to parameterize huge model i, 
we will consider out-of-core algorithm, and analyze different nu- 
merical methods to minimize ihe LSCM criterion, including multi- 
grid approaches and pre-conditioned CG. Including the stretch Ci i- 
terion directly into the LSCM criterion is also another possible fu- 
ture direction of research. 
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A PROPERTIES OF LSCMS 

The minimization problem of Section 2 has several interesting 
properties when the number p of pinned vertices in parameter spa< ;o 
is sufficient In what follows, T is assumed to be homeomorphic » 
a disc. 



A.1 Full Rank 

We first show that the matrices }A / and A have full rank whi;n 
p > 2 (p denotes the number of pinned vertices). 

For this, recall that a trjangulation that is topological^ a due 
can be incrementally constructed with only two operations (cf. Fi 5- 
ure 11): the glue operation creates one new vertex and one new 
face, and the join operation creates one new face. Thus, tnercme i- 
tal construction creates at most as much vertices as faces. Since t ie 
simplest triangulation (one triangle) has one face and three verticn «s, 
we have that n' ^ n - 2 (where n denotes ihe number of verticil 
and n' the number of triangles, as in the rest of the paper). 
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Figure 1 1 : incremental construction of a trian^ulatioti that if topoJogicatty a dl\c. 
Left: gtue two trttmghs alono an edge. Right; join two existing vertices, creating a 
new trUmgtc 



We first show thai the rank of M. / is n - p when p > 2. First 
note that since n r > t> - 2, min (n' , n - p) «= 71 - p if p £ 
2 and the rank of Mf is at most n — p. We assume that T is 
incrementally constructed with glue and join operations and prove 
the result by induction on the size of M s • We also assume, without 
loss of generality, that the p pinned vertices are concentrated in 
the initial triangulation. Let n^n* — p be the dimensions of the 
matrix at step i. Observe that since T is a non-degenerate 
Triangulation, none of the coefficients Wy,T< is zero. 

At Step 0, the triangulation has no — p =* 1 vertices and rto ^ 
1 triangles. .My 0 * has a single column and, since 7* is a proper 
triangulation, some of its coefficients are non-zero and it has rank 

~As1nimc'lnat the property holds after step L If step i + 1 is a join, 
then the number of rows of M^p grovvs by I whilo the number of 
columns is unchanged, so the rank is m+i - p = rw - p. If step 
i + 1 is a glue, a new vertex v<+i and a new triangle T arc added. 
Let vi and v? be the other vertices of T. The new matrix M f t+i) 
is as follows: 

/ I 0 \ 

I 0 

It is now easy to see that its columns are linearly independent In- 
deed, assume there are complex numbers Xj such mar 

£A>m? +1 >=0 } (5) 
J=i 

where the raj i+1) are the column vectors of Mf** 1 *. If wc 
look at the first coordinates of the column vectors, thon Equa- 
tion 5 reduces to A i m j 0 = °» which implies that A, = 
0,i = l,...,m, since M*p has full rank. Now the equa- 
tion linking tlie last coordinate of the vectors mj l+1) reduces to 
A^Wi+i.r/VSr = 0, imprying that A„, +1 = 0. Thus the 
columns of M*} +1) are linearly independent and the matrix has 
full rank. The result is proved. 

Since Mi has rank n - p, both M} and Mf have rank n - p 
when p > 2. In turn, this implies that A has rank 2(n - p) when 
p>2. 

A-2 Single Minimum 

We now show that, when p > 2, C(U) has a unique minimum. 
First, notice mat 

Now, since the rank of the Gram matrix of A (i.e. A T A) is the 
same as the rank of A A T A has rank 2(n - p) when p £ 2. Since 



-4 T -4 is a square 2(n — p) x 2(n — p) matrix, it is thus mvcrtibU 
and the minimization problem has a unique solution (when p > 2) 

x = (A T Ar 1 A T b. 

The minimum of C(U) is zero when -4x = b, i.e. when A i; 
invertible. Since it has roll rank, thiB happens exactly when A v. 
square, i.e. when n = n — p. Using the feet that n' ^ n — 2 
this implies that p = 2. Wc conclude that the mapping U is fulh 
conformal (barring self-imersections) exactly when p = 2 and thi ■ 
triangulation T is built only with glue operations. 

A.3 Invariance by Similarity 

We now prove that if U is a solution to the minimization prob- 
lem, thon zJJ + T is also a solution, for all z € C and T »■ 
(z\ . . . , z')> z ' £ C. In other words, the problem is invariant by 
a similarity transformation. 

First note that the vector H = (1, . . . f 1) T is trivially in th»i 
kernel of M> since W\ + Wi -f W'j = 0 in each triangle. Assume 
U is a solution of the problem. Wc get: 

C(zU + T) 1 zzC(\J) + 2zT*CV, 

^ ZzC(V) 4- 2z{MT)"M\J z2C(U)< 

because T m *'H is in the kernel of M. If C(V) = 0, then 
C(z\J + T) = 0. 

A.4 Independence to Resolution 

Wc now Show that if a given mesh is 'densified", then the solution l«» 
the augmented optimization problem restricted to the vertices of th : 
initial mesh is the same. Wc prove this result when 0 single triangle 
T is Split into three triangles, but the proof generalizes easily t»> 
a more general setting. So let v be the now vertex introduced in 
triangle T, i.e. as a linear combination of vertices vi, v 3| v s : 

3 3 
v = a»Vi, ^2 ai = 1, on > 0, 

Assume also for the sake of simplicity mat none of v, v J( va, v 1 
is pinned. Call Ti (i = 1, - - - ,3) the triangle created thai does nc t 
have v< as vertex. Then it is easy to see that dr* = «>rfr. 

M/ is an n' x (7* — p) matrix. After insertion or v, the ne\ > 
matrix is (n' + 2) x (n + \ - p)* Indeed, one vertex is adder , 
augmenting the number of columns by one, and three new triangles 
replace an old one, augmenting the number of rows by two. Th 5 
structure of these matrices is as follows: 




where JV) is an (n'- 1) x (n-j?) matrix. F is 1 x3> £ is 3x3 and7> 
is 3 x 1. If the coefficients of T are denoted by fj = W^.r/v 7 ^ , 
Lhcn it is easy to observe that the coefficients of C = (Uj) ani 
P=(pi) satisfy 

^ i = ^=(eti/ J -0 1> /i), pi=-^=/i, (O 
The (n -p)xl solution to the initial problem is: 
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Consider the (n ■+■ 1 - p) x 1 solution to the augmented problem: 
Uj - (Mf'MW-VMyMtU}. (7) 

Ufiing the relations of Equation 6 and the fact that /i + = 0, 
it suffices then to observe that Uj* = {Vj 7 C/v) T is the (unique) 
solution to 7, where 

f/v = oril/i -f crzUz + otaUs, 

In other words, the leasi squares conformal parameterization is un- 
changed at the old vertices and is the baryecnter of the parameteri- 
zations of vi, v 2 and vo at the new vertex v. 

A.5 Preserving Orientations 

We now sketch the proof that least squares conformal maps preserve 
orientations, i.e. there are no triangle flips. 




Figure 12: LSCMs prcsttrw orientations. 

As a preliminary, note first that if complex numbers Wi are as- 
sociated to vertices of a triangle as in Section 2.4, wiih vertices 
ordered counterclockwise, then 



(8) 



is positive (and equal to 2dr )■ 

We again assume that the tri angulation T is incrementally con- 
structed with the glue and join operations. Denote the current iri- 
angulation by Ti. For the join operation, the result is trivial. We 
now prove the result when the current step is a glue. We use the 
notations of Figure 12. Let V and V be me images ofT and 
in parameter space. Let also (resp. Wfi be complex numbers 
attached to T (resp. T') and Xj (reap. Xj) be complex numbers 
attached to V (resp. V ). Since the local bases of two triangles of 
T sharing an edge are consistently oriented, both £r — as defined in 
Equation 8 - and 

are positive. If wc assume that the unfolding of Tt has no triangle 
flips, then we also have that £ v > O f where is defined as in 
Equation 8, replacing Wj by -Xj- 

Now, writing down the equations defining U J « (U/, U*) as 
in the previous section, we find that 

WiUi + WiU* + W^Vy = 0, (9) 

where U\,Ui^ Uv arc the parameterizations of vertices vj , va, v 
(ifi, l/a being unchanged by addition of v). Using the fact that 
W[ + WZ + Wi = 0, A*J = Ui - U» and Xi ^ l/ v - Vi. 
Equation 9 rewrites as 

WlX' x - WlX f 2 = 0. (10) 
Using Equation 10 and the definition in Equation 8, we have: 
— i(X{X2 — X^Xtyi 

= 2^Li(WiWi - wiw[) = ^£Lcr' > 0 
W{W[ K 12 2 lJ WiW'i 

Thus, V' is consistently oriented and the glue operation does not 
produce a triangle flip, proving the result. 
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